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ABSTRACT 

The massive interacting binary RY Scuti is an important representative of 
an active mass-transferring system that is changing before our eyes and which 
may be an example of the formation of a Wolf-Rayet star through tidal stripping. 
Utilizing new and previously published spectra, we present examples of how a 
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number of illustrative absorption and emission features vary during the binary 
orbit. We identify spectral features associated with each component, calculate 
a new, double-lined spectroscopic binary orbit, and find masses of 7.1 ± 1.2 M 
for the bright supergiant and 30.0 ±2.1 M for the hidden massive companion. 
Through tomographic reconstruction of the component spectra from the com- 
posite spectra, we confirm the 09.7 Ibpe spectral class of the bright supergiant 
and discover a BO. 5 I spectrum associated with the hidden massive companion; 
however, we suggest that the latter is actually the spectrum of the photosphere 
of the accretion torus immediately surrounding the massive companion. We de- 
scribe the complex nature of the mass loss flows from the system in the context 
of recent hydrodynamical models for (3 Lyr, leading us to conclude RY Scuti has 
matter leaving the system in two ways: 1) a bipolar outflow from winds generated 
by the hidden massive companion, and 2) an outflow from the bright 09.7 Ibpe 
supergiant in the region near the L2 point to fill out a large, dense circumbi- 
nary disk. This circumbinary disk (radius ~ 1 AU) may feed the surrounding 
double-toroidal nebula (radius ~ 2000 AU). 

Subject headings: binaries: close — binaries: eclipsing — circumstellar matter 
- stars: individual (RY Scuti) — stars: mass loss — stars: winds, outflows 



1. Introduction 



RY Scuti (HD 169515) is a distant (1.8 ±0.1 kpc; ISmith et al.ll2002l ) and massive eclips - 
ing binary system with an orbital period currently estimated at 11.12445 days (lKreinerll2004f ). 
Analysis of the light curve indicates that at least o ne of the components fills it s Roche lobe. 



The present configuration may be semi-detached (Cowj 



( Milano et al.lll98ll ; lDjurasevic. Eshankulova. &: Erkapic 2001 ), or one with the more massive 



ev fc Hutchingd Il976l ). over contact 



comp onent embedded in an opaque thick disk (IKing fc Jamesonlll979l : lAntokhina &: Kumsiashvili 



19991 ). The binary appears t o be in an advanced stage o f evolution, and it has ejected gas 
into a young (~130 year ok 
The « 2000 AU nebula is ; 



2001 



(Gehrz et al. 


1995 


), infrared 


v images ( 


Smith et al. 


1999. 



20021 ). Since the nebula is so young and the star system is still in an active Roche- 
lobe overflow phase, RY Scuti is a powerful laboratory for studying non-spherical mass and 
angular momentum loss in interacting binaries. 



Only one of the component stars is readily y isible in optical spectra and it has a spectral 
classification of 09.7 Ibpe var (jWalbornl Il982l ). We will refer to this component as the 
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"supergiant." The other star appears to be enshrouded since it is very difficult to detect 
in the spectrum. Most investigators agree that this second star is the more massive of the 
two, and we will refer to the second component as the "massive companion." However, 
the actual masses of the stars in this system are debatable. The radial velocity curve of 
the supergiant is reasonably well established, but the results for the massive companion 
depend critically on what spectral features one assumes are associated with that star. As 
these features are difficult to ob serve and interpret, the estimated mass range has been 
huge. For example, iPopperl (119431 ) arrived at a total sys t em m ass in excess of 100M Q . Later 
investi gators found low er values: ICowley fc Hutchingd (|1976l) estimated masses of 36 and 



46M ,|SkuTskij (119921 ) found 8 and 26M , and lSahade. West, fc SkuPskvl (120021 ) estimated 



9 and 36M (for the supergiant and massive companion, respectively, in each case ) . There are 



a num ber of important photometric studies (e.g., ranging from the discover y b y iGaposch kin 



1981 



1937 through photoelectric investigations bv lGiuricin fc Mardirossian 



1981 



and Milano et al 



and up to the most recent multi-color work by iDjurasevic et al.ll200ll ). however the 
results differ with regard to the assumed binary configuration and depend sensitively on the 
mass ratio adopted from spectroscopy. 

This unique system is representative of the sh ort-lived, active mass transfer stage in the 
evolution of massive binaries. Theoretical models (jPetrovic. Langer. fc van der Huchtll2005l ) 
indicate most of the mass transfer occurs during a brief (~ 10 4 yr), and thus rare, phase in 
which the mass donor transfers most of its mass to the mass gainer star. Mass transfer will 
shrink the orbital dimensions until the mass ratio is reversed (and the gainer becomes the 
more massive star), and then the system will enter a slower (and longer lived) mass transfer 
phase as the binary expands. The massive binaries that are just emerging from t he rapid 
mass transfer phase probably belong to the observed class of W Serpentis binaries (ITarasov 
20001 ) . Only the mass donor star is visible in the spectra of these binaries and the more 
massive gainer star is hidden in a thick accretion disk (one source of emission lines in the 
spectra). The mass transfer process is complex and lea ky, and a signific ant fraction of the 
mass loss leaves the system completely (as described by iHarmaned 120021 for the best known 
object of the class, (3 Lyr). The mass donor may eventually lose its entire hydrogen envelope 
and emerge as a Wolf-Ra yet star. Therefore a system like RY Scuti may be the progenitor of 



a WR+O binary system (IGiuricin fc Mardirossian! 1 1 98 ll ; lAntokhina &: Cherepashchuklll988 
Smith et al.ll2002h . 



In a prior paper (ISmith et al.ll2002l ). several of us presented a detailed study of the 
spectral features formed in the surrounding double-toroidal nebula through an examination 
of a set of high dispersion spectra obtained with the ESO FEROS spectrograph. Here we use 
the same set of spectra supplemented by additional optical spectra to explore the spectral 
features associated with the central binary and its immediate circumstellar environment. 
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Our primary goal is to determine how the binary ejects the gas that ultimately flows into 
the dense outer double-toroidal nebula. We describe our observations and data reduction 
methods in §2 and we present in §3 examples of the orbital phase-related spectral variations 
we observed. The radial velocity curve and a new orbital solution for the supergiant mass 
donor are described in §4. Then in §5 we review the spectral clues about the nature of 
the enshrouded mass gainer, and we present a preliminary radial velocity curve and orbital 
solution for the massive companion. Both radial velocity curves are used in §6 to make 
a Doppler tomographic reconstruction of the optical spectra of the individual components. 
We describe in §7 a model for the mass outflows in RY Scuti that is based upon recent 
hydrodynamical simulations and that explains many of the observed spectral variations. 
Our results are summarized in 58. 



2. Observations and Reductions 

Our analysis is based upon spectra collected from three telescopes. The highest dis- 
persion spectra (17 in number) were obtained in 1999 with the Fiber-fed Extended Range 
Optical Spectrograph (FEROS) mounted on the 1.52 m teles cope of the European Southern 



Observatory (ESO) at La Silla, Chile fsee lSmith et alJl2002h . Also in 1999, we obtained 40 



moderate dispersion spectra (in the red region surrounding Ha) using the Kitt Peak National 
Observatory coude feed 0.9 m telescope. Finally, in 2004 we used the CTIO 1.5 m telescope 
and Cassegrain spectrograph to obtain ten blue spectra of moderate resolution covering one 
orbital period. Table [1] contains run number, dates, spectral coverage, spectral resolving 
power, number of spectra, telescope, spectrograph grating, and CCD detector used in each 
case. Exposure times were generally limited to 30 minutes or less. Each set of observa- 
tions was accompanied by numerous bias, flat-field, and ThAr comparison lamp calibration 
frames. Furthermore, we obtained multiple spectra each night of the rapidly rotating A-type 
star ( Aql for removal of atmospheric water vapor and O2 bands in the red spectra made 
with the KPNO coude feed. 

The spectra from all the telescopes were extracted and calibrated using standard routines 
111 irafS All the spectra were rectified to a unit continuum by fitting line-free regions. We 
removed the atmospheric lines from the red coude feed spectra by creating a library of ( Aql 
spectra from each run, removing the broad stellar features from these, and then dividing each 



5 IRAF is distributed by the National Optical Astronomical Observatory, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. (AURA), under cooperative agreement with the 
National Science Foundation. 
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target spectrum by the modified atmospheric spectrum that most closely matched the target 
spectrum in a selected region dominated by atmospheric absorptions. In a few cases this 
resulted in the introduction of spectral discontinuities near the atmospheric telluric lines, and 
these were excised by linear interpolation. We did not attempt any removal of atmospheric 
lines for the FEROS spectra, but some problem sections in the echelle-overlap regions were 
excised via linear interpolation. The spectra were then transformed to a common heliocentric 
wavelength grid for each of the FEROS, KPNO coude feed, and CTIO 1.5 m runs. We show 
several examples of the final spectra in the next section. 



3. Spectral Variations with Orbital Phase 



The optical spectrum of RY Scuti is very complex and time variable, and in this section 
we describe the appearance and orbital phase-related variations of several representative line 
features. A depict ion of the time-aver age of all the FEROS spectra in the range from 3600 to 
9200 A appears in lSmith et al.l (120021 ) (their Fig. 16). This illustration is ideal for identifying 
the many sharp, quadruple-peaked, emission lines formed in the surrounding double-toroidal 
nebula, but photospheric lines appear broad and shallow as the average was made over the 
full range of orbital Doppler shifts. A large fraction of the circumstellar emission lines are 
species that do not have a stellar absorption counterpart (such as [Fe III] and many other 
forbidden lines), but there are several emission lines that are superimposed upon important 
stellar spectral features (for example, in the H Balmer and He I lines). Since the nebular 



featu r es are analyzed in great detail elsewhere (IGolovatyi fe Skurskiilll992l ; ISkul'skii fe West 



19931 ; ISmith et al.l l2002f ) . we will not discuss them here. We will focus on those spectral 
features formed in the photospheres of the stellar components and in the rapidly moving 
gas immediately surrounding them. Here we introduce the different kinds of patterns of 
variability observed, and in the following sections we analyze in detail the Doppler velocity 
shifts related to features associated with the supergiant (§4) and its massive companion (§5). 

We begin with some examples from the high resolution FEROS spectra. The Si IV A4088 
absorption feature is one of only a small number of spectral lines formed in the photosphere 
of the supergiant that is not filled in or affected by nebular emission. We show the orbital 
phase variations of this line (and the nearby N III A4097 and H<5 lines) in Figure [1] as a 
function of heliocentric radial ve l ocity f or Si IV A4088. In this and the next figures we adopt 
the orbital period from iKreinerl (120041 ) of P = 11.12445 d and the epoch of phase zero as 
the supergiant superior conjunction at T$c = HJD 2,451,396.71 (derived in §4). The upper 
portion of Figure [1] shows the Si IV profiles with their continua aligned with the phase of 
observation (increasing downwards) while the lower portion shows the spectra as a gray-scale 
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image interpolated in phase and velocity. Features moving with the radial velocity curve of 
the supergiant will have a characteristic "S" shape in this image. The break in the continuity 
of the "S" curve near phase = 0.35 is due to the unfortunate gap in our phase coverage 
near there and to the simplicity of the interpolation scheme. The Si IV A4088 feature appears 
to be useful for the radial velocity measurement of the supergiant, although we find that 
the depth of the line varies with phase, weakening at = 0.0 and strengthening at = 0.5. 
There is no obvious evidence of a reverse "S" feature that would correspond to the motion 
of the massive companion (although a weak feature is present; §6). 

The orbital variations in two helium lines, the He I A4387 singlet and the He I A4471 
triplet, are illustrated in Figures [2] and [31 respecti vely. Both feature s show sharp nebular 
emission peaks from the surrounding double-torus (ISmith et al.ll2002l ) superimposed on the 
stellar absorption feature. Both lines show an absorption component that follows the radial 
velocity curve of the supergiant. However, the He I A4387 feature strengthens at both 
conjunctions and almost disappears at phase = 0.25. It also shows a blueshifted absorption 
feature that appears at = 0.5 and lasts for nearly half the orbit. Other He I singlets 
(AA4009, 4143, 4921) exhibit the same behavior. On the other hand, the He I A4471 triplet 
displays a much stronger blueshifted component between = 0.5 — 0.8 that develops into a 
very sharp and blueshifted (~ —200 km s -1 ) absorption line that lasts for greater than half 
the orbit. The other He I triplets (AA4026, 4713) also show these same features. 

Figure S] shows the orbital variations in the triplet Si III AA4552, 4567, 4574 in the veloc- 
ity frame of Si III A4552. The supergiant component appears to be present and undergoes 
the same kind of strengthening at conjunctions seen in the He I lines. However, each of the 
triplet members also show a broad shallow feature that moves in the manner expected for the 
massive companion. T hus, our observa tions confirm the detection of this second component 
that was discovered by ISkul'skiil (119921 ). 



There are two very broad emission line s present that are im portant to the study of this 
binary system, Ha A6563 and He II A4686 (jSahade et al.ll2002l ). The Ha feature consists of 
a broad emission feature that is spatially coincident with th e central binary p lus a strong 
but narrow component formed in the double-toroidal nebula ( Smith et al.ll2002l ). Therefore, 
in order to isolate the emission component near the binary, we had to remove the nebular 
components of Ha and the nearby [N II] AA6548, 6583 lines. This removal process was done 
by scaling the [N II] A6583 line in each spectrum to the appr opriate size ( of eith er Ha or 
[N II] A6548) using the equivalent widths of these lines from ISmith et al.l (120021 ). shifting 
the rescaled line to the location of the line to be removed, and then subtracting it from the 
spectrum. This process was done interactively and included small adjustments in the scaling 
and shifting parameters to optimize the subtraction. The resulting subtracted Ha profiles 
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based upon the large set of KPNO coude feed spectra are illustrated as a function of orbital 
phase in the left panel of Figure [5j We also show the one spatially resolved HST spectrum 
of the central binary from the work of ISmith et al. (which has no nebular emission 



present) that verifies that our line subtraction technique creates difference profiles with the 
appropriate shape. The emission strength appears to be much stronger at the conjunction 
phases, but this is due mainly to the drop in the continuum flux at those eclipse phases and 
our normalization of the emission strength to this varying continuum level. The right panel 
in Figure [5] shows a representation of the Ha profiles relative to a constant flux continuum 
we made by rescaling the emission flux by a factor 

///(quadrature) = KT ' 4 ^- 9 - 03 ) 

where V is the l/-band magnitud e at the orbital phase of observation (found by interpola- 



tion in the light curve data from iDjurasevic et al.ll200ll ) and V = 9.03 corresponds to the 



maximum brightness of the system at quadrature phases. The Ha emission feature appears 
broad (spanning over 1000 km s~l) and approximately constant in strength and position in 
this version. There is often a weak, blueshifted absorption component present (with a radial 
velocity of ~ —150 km s _1 ) that gives the profile a P Cygni appearance. Similar results were 
seen in the Ha difference profiles formed from the smaller set of FEROS spectra. 

Finally we show in Figure Othe orbital phase variations in the weak He II A4686 emission 
line observed in the FEROS spectra that we smoothed with a Gaussian of FWHM = 45 
km s _1 to improve the otherwise noisy appearance. This feature is found only in very hot 
plasmas, and there is no corresponding nebular feature in this case. The He II A4686 emission 
is almost as broad as the net Ha emission, but is significantly weaker. The wings of the lin e 



appear to exhibit a slight anti-phase motion that was documented by ISahade et al.l (120021 ) . 



4. Radial Velocity Curve of the Supergiant 



The radial velocity shifts of the supergiant are readily apparent in many lines, and given 
the quality and number of the new spectra we decided to measure the radial velocities and 
reassess the orbital elements. We measured relative radial velocities by cross- correlating 
each spectrum with a template spectrum. For the FEROS and CTIO spectra, this template 
spectrum was gen erated from the n o n-LT E, line-blanketed model atmosphere and synthetic 
spectra grid from lLanz fc Hubenyl (120031) using = 30 kK and log g = 3.0, which are 
appropriate values for an 09.7 lb star (IHerrero et al.lll995l ). This template was rotationally 
broadened (§5) and also smoothed to the instrumental resolution of the FEROS or CTIO 
spectra. First, we removed certain interstellar features by forming an average interstellar 
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spectrum from the mean of the entire set and then dividing each spectrum by the average 
interstellar spectrum. Because the spectrum of RY Scuti contains so many stationary nebular 
lines and possible features from the massive companion, we restricted the wavelength range 
for the cross-correlation to regions surrounding a set of absorption lines that appeared to be 
free of line blending and that were clearly visible throughout the orbit. The main features 
in the selected regions are summarized in Table [21 

For the KPNO coude feed spectra, we selected one spectrum from the set of 40 to 
serve as the template since the Lanz & Hubeny models did not match the emission lines 
we measured (Table [2]). This spectrum (made on HJD 2451425.784) has relatively high 
signal-to-noise, exhibits well-defined spectral features, and was made near quadrature when 
the supergiant Doppler shifts are large and the features reasonably well separated from any 
component from the massive companion. Once again, we first removed interstellar features 
then selected regions free from nebular lines for the cross-correlation. Once we performed 
the cross-correlation, we determined the absolute radial velocity of the template spectrum 
by fitting the core of the emission lines (N II A6610 and Si IV A6701) with a parabola and 
determined the shift from the rest wavelengths. We added the radial velocity of the template 
to all the relative velocities to place them on an absolute scale. 

Since the velocities from the different runs are based upon measurements of different 
lines in the template and individual spectra, we anticipated that there might be systematic 
differences in the zero-points of each set. We began by making independent circular orbital 
fits for each set, and indeed we found the systemic velocity 7 for the FEROS set was offset 
by —20.6 km s -1 from the resulting systemic velocities for the KPN O and CTIO sets. Since 



the latter were closer to the nebular systemic velocity (20 ±3 km s -1 ; ISmith et al.ll2002l ) and 
based upon more observations, we arbitrarily adjusted the FEROS measurements by adding 
20.6 km s _1 to bring them into consistency with the other measurements. Note that this 
decision adds an additional uncertainty to the real error in the final determination of the 
systemic velocity but otherwise has no effect on the other orbital elements. Our final radial 
velocities are collected in Table [3] that lists the date of observation, orbital phase (see below), 
radial velocity, calculated error in radial velocity, observed minus calculated residual, and 
the telescope where the spectrum was made (Tabled]). The quoted radial velocity error for 
the blue spectra is the standard deviation of the cross- correlation measurements from the 
three line regions. We measured only two line regions in the red spectra, so the velocity 
error is estimated as the larger of the difference between the two measurements or the mean 
value of \Vi — V^+i|/v2 from closely spaced pairs of observations. 

We determined revised orbital e lements using the the nonlinear, least-squares, fitting 



code of IMorbey Brosterhusl (119741 ) . Since the errors associated with each run are dif- 
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ferent, we assigned each measurement a weight proportional to the inverse squar e of the 
measu rement error. The orbital period was fixed at P = 11.12445 days as found by iKreiner 
(120041 ) from contemporary eclipse timing observations over a long time base. We made both 
eccentric and circular orbital solutions, but we think the eccentric so lution, while formally 
statistically signi ficant according to the test of iLucy fc Sweeney! (1197ll ). is probably spurious. 



Harmaned (119871 ) found that gas streams and circumstellar matter (both present in RY Scuti) 



can distort spectroscopic features in mass-transferring binaries. Such disto rtions can lead 
to skewed radial velocity results and thus artificial eccentricities (jLucyl 120051 ). A zero eccen- 
tricity is consistent with predictions for Roche lobe overflow systems where the tidal effects 
are expected to circularize the orbit and synchronize the rotational and orbital periods. The 
final orbital elements from b oth the circular and ecce ntric solu t ions a re presented in Table H] 
together with solutions from lSahade et al.l (120021 ) and lSkul'skiil (Il992l ). The epoch T sc refers 
to the time of supergiant superior conjunction (close to the time of photometric minimum 
light) while T is the epoch of periastron in the eccentric solution. The radial velocities and 
orbital velocity curves are plotted in Figure [7J Our results are in reasonable agreement with 
earlier determination s of the elements, w ith the possible exception of the larger non-zero 
eccentricity found by ISahade et al.l (120021 ) . We suspect the differ ence is due to the lack of 
orbital phase coverage in the observations of ISahade et al.l (120021 ) in the ascending portion 
of the velocity curve. 



5. Nature of the Massive Companion 

Very little is known about the massive companion because it is difficult to find its 
associated features in the spectrum (probably because the massive companion is enshrouded 
in a thick accretion disk; §7). We can make an approximate estimate for the expected 
semiamplitude of motion for the massive companion using a geometrical argument to find 
the mass ratio, q = Mmc/Msg- Let us assume the supergiant is filling its Roche lobe 
and is synchronously rotating (reasonable in t he case of RY Scuti w here we see so much 



spectroscopic evidence of active mass transfer). iGies &: Bolton! (119861 ) showed that for such 
Roche-filling stars, the projected rotational velocity V smi provides a measure of the stellar 
radius while the semiamplitude K is related to the mass ratio and semimajor axis a, and 
the ratio of these quantities is a function that solely depends on the mass ratio, 

VLRsg sin i 



V sin i/K 



M 



l+ MsG_\Rsa ={l + Q)m 
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where Q is the orbital angular rotation speed , Q = 1/q = Msg/Mmc, and f(Q) is the 
fractional Roche filling radius (lEggletonl Il983l ). Thus, a measurement of the supergiant's 
projected rotational velocity leads directly to an estimate of the mass ratio. We show below 
that this geometrical mass ratio estimate is consistent with the mass ratio derived from 
orbital velocity measurements for both components. 

We determined V sin i for the supergiant through a comparison of the photospheric line 
width with models for a grid of projected rotational velocity. For this purpose, we selected the 
Si IV A4088 line since it is free of nebular emission and line blending with other features and 
its shape is dominated by rotational broadening. We measured the FWHM of Si IV A4088 
in three of the FEROS spectra obtained near supergiant maximum velocity (to avoid the 
unusual strengthening seen at conjunctions and any line blending problems with absorption 
from the massive co mpanion). We then u sed the model profiles from the non-LTE, line 
blanketed models of lLanz fc Hubenyl (120031 ) that were rotationally broadened by a s imple 
convo lution of the zero-rotation model profiles with a rotational broadening function (IGrav 
19921 ) using a linear limb-darkening coefficient from the tables of IWade fc Rucinskil (119851 ). 
From the Lanz & Hubeny grid, we selected a model profile for T e s = 30 kK and logg = 3.0 
and adopted a linear limb darkening coefficient at this wavelength of e = 0.38. Finally, 
the resulting models were convolved with an instrumental broadening function to match 
the resolution of the FEROS spectra. The projected rotational velocity for the supergiant 
derived from the resulting (V sin i, FWHM) relation is Vsini = 80 ± 5 km s _1 . The error 
represents the standard deviation of V sin i as derived from the three measurements of line 
width and does not account for any systematic errors associated with the choice of model 
atmosphere parameters. 

We caution that this is actually an upper limit for Vsini since the lin e may also be 
broadened by macroturbulence in the stellar atmosphere (IRyans et al.l 120021). If we includ e 
an estimate of V macro = 30 km s _1 (a mid-range value for B-supergiants; IRyans et al.ll2002l ). 



then we find V eq sini = 74l 13 km s~ x (where the error range accounts for the derived error 
in Vsini and a range in macroturbulent velocity from 10 to 50 km s _1 ). Thus, using the 
V sin i/K method described above, we find an estimate of the mass ratio, q = 4.9+^. We can 
use this mass ratio estimate to offer some guidance about the expected orbital semiamplitude 
of the massive companion, K MC = K SG jq = 31 — 61 km s _1 . 

Previous researchers JPoppeiil943l ICowlev & Hutchingslfllml ISkuTskiJEiiil ISahade et al" 



20021 ) have associated a variety of spectral features with the massive companion, and we in- 
spected the orbital variations of all these proposed lines in the FEROS spectra (examples 
shown in Fig. [TJ S]) - Many of these lines have complicating factors related to nebular emis- 
sion and blending with other features. For example, a number of the He I lines appear to 
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show a blueshifted absorption feature in the phase range <p = 0.6 — 0.9 where we would expect 
to find the absorption component from the massive companion, but these blueshifted fea- 
tures first appear at (f) = 0.5 when any gas leaving the supergiant in the L2 region would also 
appear blueshifted (§7). Consequently, an identification of these features with the massive 
companion is ambiguous. 

Only three of the proposed spectral features for the massive companion clearly fit the 
mass ratio argument described above and display the expected anti-phase velocity shifts: 
Ha, He II A4686 and Si III A4552. Here we present radial velocity measurements for each 
of these. The Ha emission profile is not Gaussian in shape (see Fig. [5]) and therefore could 
not be measured as such. Instead we determined the radial velocity of the line wings (which 
represent the fastest-moving gas and are un affected by nebular emission) b a sed u pon a 



bisector position found using the method of IShafter. Szkody. fc Thorstensenl (119861 ). We 



sampled the line wings using oppositely signed Gaussian functions and determined the mid- 
point position between the wings by cross-correlating these Gaussians with the profile. We 
used Gaussian functions with FWHM = 137 km s _1 at sample positions in the wings of ±382 
km s _1 for spectra from both the coude feed and the FEROS instrument. An inspection of 
radial velocity measurements closely spaced in time indicates an average error of 1.7 km s _1 . 
Our results are given in column 3 of Table and the values are plotted as a function of 
orbital phase in Figure [71 

The He II A4686 emission profile is also broad (FWHM « 850 km s _1 ) and not Gaussian- 
shaped. Therefore, we used spectral templates and cross- correlation functions to measure 
radial velocities. We formed templates from the averages of the entire run for both the 
FEROS and CTIO data sets since the He II A4686 emission strength may have changed 
between 1999 (FEROS) and 2004 (CTIO). We then transformed the relative velocities to an 
absolute scale by adding the radial velocity derived by fitting each template with a broad 
Gaussian. Our results are given in column 4 of Table [6] and are plotted in Figure [7J Based 
on visual estimates of the goodness of the fit of this very broad line, we estimate errors in 
the velocities are ~ 50 km s _1 . 

We also measured equivalent widths for both these emission lines by a direct numerical 
integration (over the range 6536.3 — 6595.0 A for Ha and 4675.3 — 4694.5 A for He II A4686). 
Our results for Ha are given in column 4 of Table [5] and for He II are given in column 5 of 
Table [6l Errors in Ha measurements are ~0.5 A based upon examining results closely spaced 
in time. The He II equivalent width errors are 0.02 A for the FEROS data but are about 
three times worse for the CTIO data (due to the lower resolution of the CTIO spectra). The 
orbital variations of these equivalent widths are discussed below (§7). 



The absorption line Si III A4552 displays narrow and broad components that move with 
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the expected orbital Doppler shifts of the supergiant and massive component, respectively 
(Fig.HJ). In order to avoid any line blending problems between these components, we selected 
a subset of FEROS and CTIO spectra observed near the quadrature phases where the the 
broad component could be measured unambiguously, and we made Gaussian fits to determine 
the radial velocities. Using a visual estimate of goodness of fit to this broad, shallow line, 
we estimate errors in these velocities are ~ 30 km s _1 . Our results are given in column 3 of 
Table [6] and are plotted in Figure [7J 

The plots of the radial velocities in Figure [7J show that all three features exhibit the 
anti-phase motion expected for the massive companion. We made circular fits of each set 
by fixing the orbital period and epoch from the solution for the supergiant (Table H]) and 
then solving for the semiamplitude K MC and systemic velocity Jmc- The results are listed 
in Table [7J The semiamplitudes for each feature are within the expected range, however, the 
values for the two emission lines are smaller than that for the Si III A4552 absorption line. 
Furthermore, the systemic velocities for the two emission lines are larger than that found 
for the supergiant and the nebular emission. We suspect this is due to the P Cygni shape 
of both features that will bias an emission line velocity measurement to a larger value. We 
argue below (§7) that these emission features are probably formed in gas flows and are thus 
probably less representative of the orbital motion of the massive companion. 

Thus, if we adopt the Si III A4552 velocity fits as the most representative of the massive 
companion , we can make a double-lined solution of the spectroscopic masses, Msin 3 i, using 



eq. 2.52 in lffilditchl (1200 lh . 



M lj2 sin 3 i = (1.0361 x 10~ 7 )(1 - e 2 ) 3 / 2 ^ + K 2 ) 2 K 2>1 P M , 

where K is the semiamplitude in km s _1 and P is the period in days. Furthermore, if we 
set the orbital incl ination equal to that for the double-toroidal nebula, then i = 75?6 ± 1?7 



(Smit 



1981 



et al.lll999h (in subs tantial agreement with the light curve solutions; iMilano et al.l 



Djurasevic et al.ll200ll ). and we can estimate the component masses directly. We find 



the masses of the components are Msg = 7.1 ± 1.2M ffi and Mmc = 30.0 ± 2.1M M . The 



former is very low for a normal O-supergiant ([Martins. Schaerer. fc Hillier 2005 ). but we 
should bear in mind this star is in the process of losing its outer envelope to the massive 
companion and therefore its properties will be very different from those of a single star of 
comparable temperature and gravity. The mass of the massive companion is like that of 
a hot O-dwarf or an evolved B-supergiant and we will explore these and other po ssibilities 
i n the next sections. Our mass results are quit e similar to those fi rst obtained by ISkuFskii 
(119921 ) but are lower than those determined by ISahade et al.l (120021 ) (given in Table [4J . 
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6. Tomographic Reconstruction of Spectra 



Once the orbital solution fo r RY Scuti was found, we used a tomographic reconstruction 
technique (jBagnuolo et al.lll994f ) with the FEROS data to separate the individual spectra of 
the components. Tomographic reconstruction is an iterative scheme that uses the combined 
spectra and their associated radial velocities to determine the appearance of each star's spec- 
trum. RY Scuti presents a special difficulty due to the statio nary sharp nebular features. 
Therefore, before reconstruction, all nebular features listed by ISmith et al.l (120021 ) were ex- 
cised via linear interpolation. We also removed the interstellar lines from each spectrum 
prior to reconstruction to avoid creating spurious reconstructed features in their vicinity. 
The reconstruction was based upon a subset of seven FEROS spectra that were obtained 
near the velocity extrema at the quadrature phases so that we might avoid introducing ar- 
tifacts due to the line strengthening at the conjunctions and any eclipse effects. Figures [8] 
and [9] show plots of the reconstructed spectra in two different regions along wi th identifica- 



tions of the principal lines. Also plotted are two comparison spectra from the IValdes et al. 



( 120041 ) Indo-U.S. Library of Coude Feed Stellar Spectra that have a lower resolving power 
(R « 3600) than that of the FEROS spectra (R « 48000). 

1982h and we see the 



The supergiant has a classification of 09.7 Ibpe var (Walborn 



spectrum of the normal single star, HD 188209 (09.5 lab; IWalbornlll976l ). provides a rea- 
sonably good match. The He II AA4199, 4541, 4686 features appear to be slightly weaker 
in the RY Scuti supergiant, which is consistent with its subtype difference in spectral type. 
We used the ratios of the equivalent widths of several weaker He I lines for the RY Scuti 
supergiant and HD 188209 spectra to set the monochromatic flux ratio for the RY Scuti 
binary components, and we find the supergiant contributes ~ 57 ± 13% of the total flux. 
The error is based on scatter between the results from different spectral features. The re- 
constructed spectra are plotted in Figures [8] and M normalized to their respective continuum 
fluxes, so individual lines appear deeper than in the composite spectra where each is diluted 
by the flux of the other star. There are several obvious differences between the supergiant 
component and single star spectrum in some of the stronger lines such as He A3970 and 
US A4101 but we do not ascribe any particular significance to these as such features are 
clearly distorted by emission from circumstellar gas. On the other hand, there does appear 
to be a significant weakness in the C lines and an enhancement in the N lines in the spec- 
trum of the RY Scuti supergiant compared to that of HD 188209. For example, the C III 
AA4070,4647 lines that are present in the spectrum of HD 188209 are conspicuously absent 
in the RY Scuti supergiant's spectrum. This suggests the atmosp here of the superg iant is 
enhanced with CNO-processed gas (as is the surrounding nebula; ISmith et al.ll2002l ). Per- 
haps this is not surprising given that a large portion of the supergiant's mass must have 
already been lost to reveal gas from deeper layers closer to the source of core H-fusion. This 



-14- 



result suggests that the bright 09.7 supergiant is the ultimate source of the processed gas 
that now forms the outer double-toroidal nebula. 



Skul'skiil (Il992l ) suggested the spectral features moving with Doppler shifts of the mas- 



sive companion have the appearance of a B2 star. We find the lines in the reconstruction 
for the orbital sh ifts of the massive companion mo re closely resemble that of the BO. 5 la 
star, HD 185859 (Morgan. Code, fc Whitfordlll955l ). However, there are several features in 
the reconstruction that make us caution against assuming this spectrum forms in the pho- 
tosphere of a supergiant star. First, the luminosity sensitive lines of Si IV AA4088,4116 and 
Si III AA4552, 4567, 4574 appear very strong in the reconstructed spectrum relative to those 
in HD 185859 which would imply the lines fo rm in a very low density plas ma (like those in 
the photospheres of the most luminous stars; IWalborn fc Fitzpatricklll990l ). Secondly, if we 
suppose that the gravity of a BO. 5 la star is logg = 3.0 (see the case of the like star k Ori; 
McErlean. Lennon. fc Duftonl 119991 ). then using our derived mass, the massive companion 
star would have a radius of ~ 29R Q , which is 83% of its Roche radius. A radius this big 
would have two observational consequences: first, the flux from this component would be 
twice the flux from the bright 09.7 supergiant, and second, a star this large would cause the 
spectral lines of the bright 09.7 supergiant to weaken or disappear during primary eclipse 
(<ft = 0.0). Neither of these two predictions are seen: the massive companion flux is less than 
that of the supergiant and the supergiant's absorption lines do not weaken at primary eclipse 
(Fig. []]- HI). Also, we estimate the projected rotational velocity of the massive companion 
spectrum by measuring the widths of the Si III A4552 feature in five FEROS spectra where 
the feature was well separated from any component of the bright supergiant. We used these 
widths to obtain V sin i = 194 ± 18 km s" 1 in the same way as described above for the pri- 
mary (§5) by comparing them with rotationally broadened model spectra for T e g = 27500 K, 
log q = 3.0, and a l imb darkening coefficient e = 0.34. This result is unreasonably large since 



Ryans et al.l (120021 ) find that no early B-type supergiant has a V sin i greater than 60 km s -1 . 
Therefore, we suggest these characteristics are best explained if we place the origin of the 
spectrum in the photosphere of a thick accretion disk surrounding the massive companion 
rather than in the photosphere of the star itself and we explore this idea further in the next 
section. Note that the very broad He II A4686 emission is assigned to the massive companion 
in the reconstruction, but we argue below that its site of origin may not be exactly co-spatial 
with that of the absorption line spectrum. 
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Mass Outflows in RY Scuti 



Many of the spectral features discussed above originate in gas outflows from the binary, 
and in this section we explore their sites of origin with reference to recent hydrodynamical 
simulations of gas flows in Roche lobe overflow binaries. The luminous supergiant in RY Scuti 
can potentially be losi ng mass by a radi a tively driven wind and/or tidal streams along the 
axis joining the stars. IFriend fc Castor! (119821 ) show that a continuum of states can exist 
between these idealized cases. When a luminous star has a radius that is significantly smaller 
than the critical Roche lobe, the wind will be fast and spherically symmetric (as for single 
stars). As the ratio of stellar to Roche radius increases, the star becomes tidally distorted 
and its wind becomes asymmetric with a higher mass loss rate and slower outflow along the 
axis joining the stars (a "focused wind"). Finally, once the star essentially attains a size 
filling the Roche surface, most of the mass loss will occur in slow outflows directed along the 
axis joining the stars. All of the photometric light curve studies confirm the supergiant is 
very tidally distorted and must therefore be close to Roche filling. Thus, we will explore the 
consequences of tidally formed outflows in reviewing the observational clues about the mass 
loss. 

The discussion begins with the dimensions of the system, and in Figure [10] we present 
a cartoon illustration of the cross section of the system as viewed from above the orbital 
plane. The direction of the observer at different orbital phases is indicated by numbers in 
the periphery of the diagram. As before, we assume the supergiant fills its criti cal Roche 
surfa ce and thus has an equivalent volume radius of approximately Rsg = 18i? (lEggleton 
19831 ) . The supergiant is depicted on the left side of Figure [KB If we assume the massive 
companion is a main sequence star, then for its derived mass we expect it is an 6.5 V 
star with a polar radius of approximately 9.QR & (see Table 4 in iMartins et al.l 120051 ) . The 
equatorial radius could be significantly larger if the star is rapidly rotating, but we assume 
it is a spherical star in its depiction on the right hand side of Figure [KB Either way, the 
radius of the massive companion is probably significantly less than its critical Roche radius, 
and therefore we begin our analysis assuming that the system is semi-detached. 

The kinds of gas outflows expected in s ystems like RY Scuti are explored in a series 
of recent papers bv iNazarenko fc Glazunov^l (120031 . l2006al lbl) . They present hydrodynamical 
simulations in two and three-dimension s to model the case of the interacting binary (3 Lyr 
(jHarmanec et al.lll996l ; lHarmaned 120021 ) . RY Scuti and f3 Lyr share many of the properties 
common to the W Serpentis class, and we can obtain considerable insight about the gas flows 
in RY Scuti from these papers (although we caution that the stars in RY Scuti are both hot- 
ter and more massive than those in (3 Lyr). These numerical models assume the mass donor 
is a low gravity object that fills its Roche surface, and the authors use the models to follow 
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the time evolution of the gas flows until they reach a semi-static configuration. There are 
several features of these models that are relevant to the case at hand. First, Nazarenko & 
Glazunova find that the donor loses mass through both the inner LI region and the outer L2 
region (facing away from the companion). Because the gas is so loosely bound to the donor, 
these gas streams are not compressed into a classical stream but instead fan out across a 
significant portion of the stellar surface facing towards and away from the companion. The 
approximate ranges of these outflows are indicated by the dashed lines to the left and right of 
the supergiant in Figure [TD] and are based upon simple ballistic traject ories. The thicker lin e 
to the right of the supergiant indicates the classical LI gas trajectory (ILubow fc Shulll975l ). 
The models also show an accretion disk does form but that gas may be lost from the binary 
near the outer L3 point (to the right of the massive companion in Fig. [TOl) . We have simply 
illustrated the disk in Figure (TU] assuming it extends to 80% of the Roche radius of the mas- 
sive co mpanion (RR OC he ~ 35i? ). The three-dimensional models (INazarenko fc Glazunova 
2006al lbh show the inner portions of the disk attain a significant height above and below the 



orbital plane, and therefore the gas structure is referred to as a torus rather than a disk. If 
the massive companion is a main sequence star, then the height of the torus may well exceed 
the st ellar radius and block it from vi ew for all but extremely low orbital inclinations. Fi- 
nally, INazarenko fc Glazunoval (l2006bl ) consider the effect that a stellar wind from the mass 
gainer will have on the surrounding torus. They find the wind breaks out in bipolar outflows 
and that the hottest gas temperatures occur where the wind strikes the inner torus. Such a 



bipol ar outflow may explain the elongated radio emission surrounding f3 Lyr (jUmana et al. 



20001 ) . If the massive companion in RY Scuti is an 06.5 V star, then we expect it will be 



hot and luminous (perhaps exceeding the luminosity of normal dwarfs due to the addition 
of accretion luminosity) and thus it will be a source of a significant radiatively-driven wind. 

We will use the framework of gas features illustrated in Figure [10] to offer a plausible 
origin for some of the spectral features discussed in earlier sections. We begin with the 
expected outflows from the supergiant towards LI and L2. Most of the absorption lines 
associated with the supergiant appear to strengthen considerably in depth at the conjunction 
phases and this change is particularly striking in the case of the He I lines illustrated in 
Figures [2] and [3J The absorption preferentially strengthens on the blue side of the profile, 
and in the case of supergiant inferior conjunction {(f) = 0.5), remains blueshifted for a while 
longer. For triplet transitions like He I A4471 (Fig. [3]), the blueshifted component transforms 
from a broad feature into a narrow, almost nebular absorption line that lasts for more than 
half the orbit. We suggest these changes are best explained by absorption from the outflows 
from the supergiant that emerge in the LI and L2 directions. The blueshifts we observe with 
projected speeds of ~ 200 km s _1 are consistent with an outflow towards LI (at = 0.0) and 
towards L2 (at <fi = 0.5). The disappearance of the broad absorption shortly after (f> = 0.0 
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is expected since this flow terminates only a short distance from the supergiant where the 
gas strikes the accretion disk. On the other hand, the gas leaving the L2 region forms a long 
trailing plume that may occult the supergiant for a considerable time following <fi = 0.5 and 
thus may account for the presence of blueshifted absorption until <fi ~ 0.8. The appearance 
of the narrow blueshifted absorption in the He I triplets may be due to low optical depth 
gas associated with the L2 outflow that extends well above and below the orbital plane. 

The gas that enters the accretion disk from the LI outflow is probably very dense and 
extended from the orbital plane. In such a situation, the massive companion may be wholly 
obscured from our point of view and its photospheric spectrum will be invisible. However, 
we noted in §5 and §6 that there are a number of spectral features that do display radial 
velocity patterns opposite to those of the supergiant and may plausibly be associated with the 
immediate environment of the massive companion. In particular, we suggest the Si III and 
other lines found in the spectral reconstruction for orbital motion of the massive companion 
(§6) are probably formed in the extended accretion disk "photosphere" surrounding the 
massive companion. These features are very broad (indicative of the large Keplerian motions 
in the disk) and appear similar to those in early B-supergiants (suggesting a gas density that 
is relatively low compared to most stellar photospheres). Furthermore, the fact that this 
spectral component appears slightly cooler (like a BO. 5 I star) than the supergiant (09.7 Ibpe) 
is consistent with idea that the gas has cooled slightly since leaving the supergiant and is 
shielded from the flux from the presumably hot massive companion by the denser regions of 
the accretion disk. Simila r kinds of disk photospheric features have recently been detected 



in the spectrum of (3 Lyr (lAk et al.ll2007l ). 



Minimum light according to the recent ephemeris from iKreinerl (12004 ) occurs at spec- 
troscopic phase = 0.980 ± 0.008 (Table HJ). The fact that the deeper eclipse occurs when 
the supergiant is behind the massive component is consistent with the somewhat cooler disk 
photosphere we find from the line patterns in the reconstructed spectrum. 

The large mass of the massive companion suggests the star probably has a vigorous 
radiatively-driven stellar wind. Hydrodynamical models for f3 Lyr indicate the radiatively- 
driven wind of the mass gainer will break out of the accretion torus in bipo l ar flows that ar e 



the suspected source of the Ha emission in j3 Lyr (lHarmanec et al.l Il996l ; lAk et all 120071 ) . 
If bipolar wind outflows are also the source of some of the Ha emission in RY Scuti, then 
they should extend far away from the orbital plane and be minimally eclipsed at conjunction 
phases. We may test this idea by comparing the orbital phase variations in the Ha equivalent 
width (W\) with the l/-band light curve. We measure the total emission flux relative to 
a time variable continuum flux and if the emission flux is constant, then the product of 
the emission equivalent width and continuum flux will also be constant (i.e., the measured 
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equivalent width will vary inversely as the continuum flux). In Figure [TTJ we show the 
orbital phase variations in Wx(Ha) normalized to their mean at the quadrature phases. We 
see the measured W\ increases during the decreases in flux that occur surrounding the times 
of eclipse at <fr — 0.0 and 0.5. We also sho w in Figure [TT1 the i nverse continuum flux variation 



derived from the l^-band observations of iDjurasevic et al.l (120011 ) and transformed by the 



equation described in §3. The fact the broad Ha equivalent width variations follow the 
inverse light curve so closely indicates the Ha emission source is not substantially eclipsed 
by either the supergiant or the accretion disk, and therefore probably forms over a volume 
that is either very large and/or extended far from the orbital plane (we discuss this further 
below). We showed in §5 that the bisector velocities of the Ha emission wings (corresponding 
to the fastest moving gas) show approximately the same velocity pattern as expected for the 
massive companion and this anti-phase velocity pattern is seen in Figure [5J The combined 
evidence of orbital motion and the lack of eclipses suggests the emission line wings of the 
Ha feature (and probably part of the main core) form in the bipolar outflows from the disk- 
constrained wind of the massive companion. The highest projected speeds observed (~ 500 
km s _1 ) are consistent with outflow velocities normal to the orbital plane of ~ 2000 km s _1 , 
a typical value among O-dwarf winds. 

The central part of the Ha profile appears as a quasi-P Cygni profile (with an emission 
peak at redshifted velocities and a weak absorption at blueshifted velocities) that usually 
implies a gas outflow. This part of the Ha profile is generally stationary, showing neither 
the motion associated with the supergiant nor that of the massive companion. This is very 
different from the case of the wind emission of the comparable supergiant in the massive 
X-ray binary system Cygnus X-l, w here the Ha em ission is much weaker and follows the 



Doppler motions of the supergiant (IGies et al.l 120031 ). Furthermore, the wind emission in 



single OB-supergia nts is generally an order of magnitude weaker than observed in RY Scuti 



(IMorel et al.ll2004j ). which indicates that the dense outflow in RY Scuti occurs over a volume 
much larger than that in the region immediately surrounding the supergiant. The lack of 
any orbital velocity variation suggests either the source is close to the center of mass or 
the flux originates far from the center of mass where any orbital motions disappear due to 
conservation of angular momentum. The first possibility appears to be ruled out by the fact 
that the main core shows no evidence of eclipses (Fig. ITTT) . Taken together, these properties 
suggest that a large fraction of the emission flux originates in outflowing gas that surrounds 
the binary as a whole. This circumbinary disk is probably fed by the gas flows leaving 
the binary in the vicinity of the L2 and L3 regions. These broad flows may have sufficient 
vertical extent to account for the blueshifted absorption observed at the quadrature phases. 
The absence of any substantial decrease in emission flux around <fi = 0.0 (Fig. [LTJ suggests 
that the inner radius of the circumbinary disk is larger than asc + Rsg/ cosi = 131i? , i.e., 
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approximately twice the binary separation. On the other hand, since the gas density will 
decline with increasing radius and since the Ha emission depends on gas density squared, 
we expect that most of the emission comes from the inner part of the circumbinary disk. 
Thus, the Ha emitting region of the circumbinary disk probably has a characteristic radius of 
approximately 1 AU. The circumbinary disk flux contribution appears to decline as expected 



among the profiles of the higher members of the Balmer sequence (see Fig. 16 in lSmith et al. 



20021 ). We note for completeness that the broad asymmetrical shape of the Ha emission 
from the inner circumbinary disk appears som ewhat like a veloc ity-expanded version of the 
profile from the outer double-toroidal nebula (jSmith et al.ll2002l ). Since the latter is due to 
a spatial gas asymmetry in the outer expanding ring, it is possible that there may also be 
asymmetries in the mass loss into the inner circumbinary disk. 

The He II A4686 emission found in the spectrum of RY Scuti has a radial velocity curve 
that is similar to that of the massive companion (Fig. [TJ and this suggests it is formed in 
the same vicinity. Furthermore, the equivalent width variations of the He II A4686 emission 
(Fig. [TTj) indicate a drop in relative emission strength at <fi = 0.5 when the supergiant is in 
the fore ground, consistent with an origin near the massive companion. The hydrodynamical 
models (INazarenko fc Glazunovall2006bl ) suggest that the hottest gas is located at the place 
where some of the incoming gas stream through LI climbs over the accretion torus vertical 
extensions and collides with the high speed wind from the massive companion. Thus, the 
hot gas is predicted to have a velocity curve for a position slightly offset from the centroid 
of the massive companion towards the binary center of mass. This appears to be consistent 
with the He II A4686 velocity curve (Fig. [7J) that has a slightly smaller semiamplitude than 
that for the Si III lines that probably form in the disk photosphere (which is approximately 
centered on the massive companion). 

Our results indicate the mass outflow from RY Scuti has at least two components, 
a bipolar outflow of hot stellar wind gas from the massive companion and an equatorial 
outflow of cooler gas fed by the plume from the L2 region and disk leakage from the L3 
region. There is no direct evidence of the high-speed, bipolar outflow in Ha images made 
with HST (shown in ISmith et al.lll999l ). but future optical and radio emission maps made 
with much higher angular resolution may reveal these jet-like features. We suspect most of 
the cool gas outflow in the orbital plane ends up in a circumbinary disk that slowly transports 
mass and angular momentum away from the binary and may eventually be a source of gas 
for the outer double-toriodal nebula. This draining of angular momentum fr om the binar y 
may explain the observed shrinking of the orbital period and semimajor axis (jSimonlll999l ). 
It is possible the density and spatial extent of the circumbinary disk is sufficient to become 
a significant opacity source blocking the ionizing flux from the central binary. We speculate 
that this may be one of the reasons why the surrounding H II double-toroidal nebula seen in 
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HST images has a darker equatorial zone (shaded from ionizing radiation by the intervening 
circumbinary disk). 



8. Conclusions 

Our spectroscopic investigation of RY Scuti has led to a reassessment of the component 
masses, a new understanding of the properties of the spectral features related to each star, 
and a sketch of how the binary is ejecting gas from the system. Our new radial velocities 
and orbital elements for the supergiant improve upon the past work and we argue that 
the orbit is probably circular as expected for Roche-filling systems. We determined radial 
velocities for the Si III A4552 absorption line, a feature that is probably formed in the 
vicinity of massive companion, and by combining the semiamplitudes of the supergiant and 
this feature, we estimate the mass ratio is Mmc/Msg = 4.2 ± 0.7. If we adopt an orbital 
inclination from that for the surrounding double-toroidal nebula, then the stellar masses 
are Msg = 7.1 ± 1.2M and Mmc — 30.0 ± 2.1M Q . These are consistent with the idea 
that the massive companion is an O-star hidden within a dense accretion disk and that the 
mass-losing supergiant is destined to become a Wolf-Rayet star. In the distant future, this 
system may resemble the well known O + WR bi nary t 2 Vel that h as component masses of 



M(0) = 28.5 ± 1.1M Q and M(WR) = 9.0 ± 0.6 flNorth et al.l 120071 ) 



We used a Doppler tomographic reconstruction scheme to separate the component op- 
tical spectra of the stars. The line patterns of the the supergiant dominate the combined 
spectrum and a comparison with a similarly classified star indicates that the supergiant 
contributes about 57% of the visible flux. The optical spectrum of the supergiant displays 
unusually weak carbon lines, suggesting the photosphere (like the surrounding nebula) con- 
tains CNO-processed gas. The reconstructed spectrum for the massive companion has a 
superficial resemblance to that of a rapidly rotating B-supergiant, but we argue this spec- 
trum probably forms in the accretion disk photosphere surrounding the massive companion. 

We found several lines of evidence indicating mass loss from the binary. The apparent 
strengthening of blueshifted parts of many absorption lines at both conjunctions and espe- 
cially following supergiant inferior conjunction suggest that the supergiant is losing mass in 
broad streams directed both towards and away from the massive companion. The strong 
Ha feature appears as broad P Cygni line indicating outflow. The radial velocity curve of 
the Ha wings indicates the highest velocity emission originates near the massive companion, 
probably in a bipolar outflow where the stellar wind of the massive companion breaks out 
of the thick surrounding accretion disk. The core Ha emission is stationary and uneclipsed, 
and we suggest it forms mainly in an outflowing equatorial region (a circumbinary disk). 
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The other important emission line, He II A4686, also has a radial velocity curve consistent 
with an origin near the massive companion, and we suggest it forms in hot regions where 
gas from the LI stream flows over the thicker portions of the disk and strikes the fast wind 
from the massive companion. 

We argue that the gas lost from the supergiant through the outer L2 region and also 
lost from the outer parts of the accretion disk near L3 ends up in a relatively cool and dense 
circumbinary disk with a characteristic Ha emitting radius of 1 AU. We suspect this disk gas 
will con tinue to flow outwar ds and may help channel gas into the surrounding double-toroidal 
nebula. ISmith et al.l (120011 ) studied the proper motions of the outer double-torus nebula and 
they conclude the gas in this H II region was probably ejected from the central binary around 
the year 1876 ± 20. It is possible this time corresponded to an epoch with an especially high 
binary mass transfer rate and associated outflow int o the circumbinary disk. We note, for 
example, that the first light curve for the binary from I G ap oschkinl (119371 ) (covering the years 
1924 to 1934) shows a second ary minimum that is on ly 0.3 mag deep compared to the current 
day observed 0.5 mag drop (IDjurasevic et al.l 120011 ) while the primary eclipse depth is the 
same. This suggests the disk surrounding the massive companion was cooler in the early part 
of the last century, which would be consistent with the presence of a thicker, more extended 
disk resulting from a higher mass transfer rate. Finally, we note that the mass loss into 
the circumbinary disk removes angular momentum fro m the centra l binary, consistent with 
the observation that the orbital period is decreasing (jSimonlll999l ). Thus, contrary to the 
predictions of conservative mass transfer models, the components of RY Scuti are continuing 
to be drawn together even after the mass ratio reversal has occurred. We may therefore 
expect that an episode of more intense interaction and mass transfer and mass loss lies in 
the immediate future for RY Scuti, and continuing observations will be vital in documenting 
this key stage in the evolution of massive binaries. 
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Table 1. Journal of Spectroscopy 



Run Dates Range Resolving Power Number of Observatory/Telescope/ 

Number (HJD-2,450,000) (A) (A/AA) Spectra Spectral Grating/CCD 



1 1354.7- 1354.8 6431 - 6785 5440 2 KPNO/0.9m/RC181/TI5 

2 1355.7 - 1364.9 5405 - 6743 3950 32 KPNO/0.9m/RC181/F3KB 

3 1373.7 - 1394.6 3600 - 9200 48000 17 ESO /1.5m/FEROS/EEV 2 x 4K 

4 1421.8 - 1429.7 5397 - 6735 4050 6 KPNO/0.9m/RC181/F3KB 

5 3152.9 - 3164.9 4068 - 4738 2430 10 CTIO/1.5m/#47II/Loral IK 



Table 2. Supergiant Radial Velocity Line Sample 



ESO 1.5 m Si IV A4088; Si IV A4116, He I A4121; 

NIII AA4510, 4514, 4518; 
CTIO 1.5 m Si IV A4088; Si IV A4116, He I A4121; 

NIII AA4510, 4514, 4518; 
KPNO CF N II A6610 (em.); Si IV A6701 (em.) 
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Table 3. Supergiant Radial Velocity Measurements 



HJD 


Orbital 


V r 


a 


O — C 




(—2,450,000) 


Phase 


(km s L ) 


(km s L ) 


(km s 1 ) 


Telescope 


1 354 781 


0.231 


—221.6 


11.6 


14.6 


KPNO CF 


1 354 SOS 

-LiJiJ^. OUO .... 


0.233 


—218.0 


9.6 


18.5 


KPNO OF 


1 355 fi74 


0.312 


—217.2 


30.0 


2.2 


KPNO OF 


1 355 701 


0.314 


—207.2 


11.3 


10.8 


KPNO CF 


1 355 Q1 4 


0.333 


—215.4 


16.9 


—10.8 


KPNO CF 


1 355 035 


0.335 


—212.4 


27.4 


—9.3 


KPNO CF 


1 356 722 

-LiJtJU. 1 iiii .... 


0.406 


—116.8 


20.1 


10.9 


KPNO OF 


1 356 744 

-LiJtJU. 1 Ttrfc .... 


0.408 


— 122.4 


25.1 


2.6 


KPNO OF 


] 356 qin 


0.423 


— 122.0 


28.7 


—17.0 


KPNO CF 


1 357 691 

4Ul/ 1 .V^C-L • . • . 


0.493 


—30.1 


12.2 


—30.5 


KPNO CF 


1 357 71 2 

UU 1 . f J-iJ .... 


0.495 


—7.8 


10.1 


— 11.1 


KPNO CF 


1357.832 


0.506 


8.5 


13.4 


-11.7 


KPNO CF 


1357.853 


0.507 


25.2 


34.2 


2.0 


KPNO OF 


1357.912 


0.513 


21.5 


9.7 


—9.9 


KPNO CF 


1 357 Q33 

-LiJfJ I . <7<J'-J .... 


0.515 


29.7 


9.6 


—4.7 


KPNO OF 


1 358 8fi1 

_1_ i_) tj (J . O <J -L .... 


0.598 


127.5 


9.6 


—28.1 


KPNO OF 


1 35q 663 


0.670 


232.0 


26.9 


2.0 


KPNO OF 


1 35a 70Q 


0.674 


214.8 


9.8 


— 18.3 


KPNO CF 


i 35a ssq 

lUt^C.UUC .... 


0.691 


213.2 


13.3 


—30.4 


KPNO CF 


1 35q 01 1 

4 U U f7 i C 1 4 .... 


0.692 


223.3 


15.7 


—21.4 


KPNO CF 


1 3fi0 fifi4 

lUUU.UU^ .... 


0.760 


251.5 


10.7 


—8.8 


KPNO OF 


1 Sfifl 7,30 


0.766 


248.9 


17.6 


— 10.6 


KPNO CF 


1 360 SfiO 

iuuu,uu w .... 


0.778 


238.7 


25.9 


— 18.3 


KPNO CF 


1 360 881 

UUUtUUl .... 


0.780 


242.1 


27.4 


— 14.4 


KPNO CF 


1 361 »72 


0.869 


209.2 


20.6 


14.7 


KPNO CF 


1 3fil 953 

_L ij \J J_ . C '-J .... 


0.876 


199.8 


10.0 


13.2 


KPNO OF 


1 362 661 

4'JUiJ , UU 4 . . . . 


0.940 


142.8 


9.6 


39.1 


KPNO CF 


1 362 706 

lOUi. 1 LJU .... 


0.944 


120.8 


15.2 


22.9 


KPNO OF 


loOO.DO^t .... 


U.UOZ 


U.O 


y.o 


oo.o 




1363.807 .... 


0.043 


-43.5 


9.6 


11.0 


KPNO CF 


1363.839 .... 


0.046 


-61.5 


9.6 


-2.6 


KPNO CF 


1363.925 .... 


0.053 


-70.6 


9.6 


-0.1 


KPNO CF 


1363.946 .... 


0.055 


-75.4 


9.6 


-2.1 


KPNO CF 


1364.900 .... 


0.141 


-191.6 


31.9 


-10.1 


KPNO CF 


1373.693 .... 


0.931 


112.4 


12.0 


-3.3 


ESO 1.5 m 


1373.758 .... 


0.937 


103.7 


8.6 


-3.7 


ESO 1.5 m 


1374.655 .... 


0.018 


-53.7 


7.3 


-37.3 


ESO 1.5 m 


1375.709 .... 


0.113 


-161.2 


6.5 


-10.6 


ESO 1.5 m 


1379.718 .... 


0.473 


-16.1 


7.4 


14.6 


ESO 1.5 m 


1381.646 .... 


0.646 


214.1 


11.5 


4.4 


ESO 1.5 m 


1382.668 .... 


0.738 


272.0 


3.9 


11.9 


ESO 1.5 m 


1383.645 .... 


0.826 


252.3 


8.5 


19.4 


ESO 1.5 m 


1384.625 .... 


0.914 


140.1 


6.8 


0.5 


ESO 1.5 m 


1385.687 .... 


0.010 


-53.4 


10.9 


-49.9 


ESO 1.5 m 


1386.650 .... 


0.096 


-161.5 


6.9 


-31.5 


ESO 1.5 m 
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Table 3 — Continued 



HJD 


Orbital 


V r 


a 


O-C 




(-2,450,000) 


Phase 


(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 


Telescope 


1388.583 .... 


0.270 


-215.0 


3.4 


20.9 


ESO 1.5 m 


1390.664 .... 


0.457 


-49.4 


8.9 


5.8 


ESO 1.5 m 


1391.602 .... 


0.541 


78.5 


12.8 


3.2 


ESO 1.5 m 


1392.619 .... 


0.633 


194.5 


10.4 


-1.6 


ESO 1.5 m 


1393.638 .... 


0.724 


272.4 


7.1 


14.9 


ESO 1.5 m 


1394.615 .... 


0.812 


253.2 


1.4 


11.2 


ESO 1.5 m 


1421.765 .... 


0.253 


-230.9 


9.6 


6.9 


KPNO CF 


1425.784 .... 


0.614 


165.2 


9.6 


-9.9 


KPNO CF 


1426.745 .... 


0.700 


236.6 


29.3 


—12.1 


KPNO CI 


1427.748 .... 


0.790 


244.5 


27.8 


-8.2 


KPNO CF 


1428.718 .... 


0.878 


179.7 


24.1 


-5.0 


KPNO CF 


1429.713 .... 


0.967 


98.4 


21.3 


35.8 


KPNO CF 


3152.930 .... 


0.871 


195.5 


13.9 


3.0 


CTIO 1.5 m 


3153.825 .... 


0.951 


83.9 


22.5 


-3.0 


CTIO 1.5 m 


3154.927 .... 


0.050 


-118.3 


46.0 


-52.4 


CTIO 1.5 m 


3156.922 .... 


0.229 


-216.3 


15.4 


19.5 


CTIO 1.5 m 


3157.932 .... 


0.320 


-204.4 


17.0 


9.6 


CTIO 1.5 m 


3158.928 .... 


0.410 


-128.2 


2.5 


-5.9 


CTIO 1.5 m 


3161.915 .... 


0.678 


239.7 


14.6 


3.8 


CTIO 1.5 m 


3162.914 .... 


0.768 


254.1 


11.2 


-5.0 


CTIO 1.5 m 


3163.892 .... 


0.856 


206.1 


17.7 


-1.4 


CTIO 1.5 m 


3164.903 .... 


0.947 


83.5 


19.6 


-9.5 


CTIO 1.5 m 



Table 4. Orbital Elements for RY Scuti 



Element 


Circular 


Elliptical 


Sahadc ct al. (2002) 


Skul'skii (1992) 


P (days) 


11.12445 a 


11.12445 a 


11.124646 a 


11.1250 a 


T sc (HJD-2,400,000) 


51396.71 ±0.02 






44777.9 ± 1.0 


T (HJD-2,400,000) . . 




51395.68 ± 0.40 


45107.74 




e 





0.05 ±0.01 


0.16 ±0.04 





w (deg) 




57 ± 13 


345 ± 15 




K S g (km s _1 ) 


249 ± 3 


247 ±3 


273 ± 14 


235 ±3 


Kmc (km s _1 ) 


59 ±9 




71 ±9 


71 ±2 


7SG (km s- 1 ) 


11 ±2 


12 ±2 


17 ±5 


9±2 


IMC (km s _1 ) 


-20 ±9 






11±2 


r.m.s.jg (km s _1 ) ... 


19 


17 






r.m.s.jKC (km s _1 ) .. 


25 








a (Rq) 


69.9 ±2.3 




77 ±4 


68 


q-M MC /M SG 


4.2 ±0.7 




3.9 ±0.5 


3.3 ±0.1 


M SG (Mq) 


7.1 ± 1.2 




10.2 ±1.7 


8 


M M C (Mq) 


30.0 ±2.1 




39 ±4 


26 



a Fixcd 
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Table 5. Ha Wing Velocities and Equivalent Widths 



Date 


Orbital 


V r 


W x 


(HJD — 2,450,000) 


Phase 


(km s L ) 


(A) 


1354 781 


0.231 


83.3 


— 13.4 


1 354 803 


0.233 


82.9 


— 13.4 


1 355 674 


0.312 


76.8 


— 12.9 


1 355 7(11 


0.314 


73.4 


— 13.2 


1355 913 


0.333 


74.2 


— 13.2 


1 355 935 


0.335 


72.3 


— 12.7 


1 35fi 722 


0.406 


60.3 


— 14.8 


1 356 744 


0.408 


61.0 


— 14.8 


1356 910 


0.423 


59.0 


— 14.6 


1357 691 


0.493 


42.6 


—21.9 


1357 712 


0.495 


44.1 


—22.8 


1357 832 


0.506 


38.8 


—23.0 


1 357 853 


0.507 


38.0 


—23.0 


1357 912 


0.513 


35.7 


—20.5 


1 357 933 


0.515 


35.4 


— 19.6 


1 358 861 


0.598 


30.6 


— 14.5 


1 359 663 


0.670 


20.8 


— 12.5 


1 359 709 


0.674 


23.7 


— 13.4 


1 359 889 


0.691 


21.9 


— 12.1 


1359 911 


0.692 


21.8 


— 12.3 


1360 664 


0.760 


26.4 


— 13.6 


1 360 730 


0.766 


26.3 


— 14.0 


1 360 860 


0.778 


33.4 


— 13.0 


1360 881 


0.780 


33.8 


— 12.7 


1361 872 


0.869 


38.5 


— 13.8 


1361 953 


0.876 


40.2 


— 12.1 


1362 661 


0.940 


27.3 


— 18.8 


1 362 706 


0.944 


25.7 


— 19.5 




U.UOZ 


AQ ^ 


91 Q 


1363.807 


0.043 


55.2 


-21.3 


1363.839 


0.046 


56.9 


-21.1 


1363.925 


0.053 


66.0 


-18.6 


1363.946 


0.055 


70.3 


-18.0 


1364.900 


0.141 


119.1 


-13.6 


1373.693 


0.931 


24.4 


-18.4 


1373.758 


0.937 


24.4 


-19.2 


1374.655 


0.018 


39.5 


-21.5 


1375.709 


0.113 


90.6 


-15.1 


1379.718 


0.473 


41.7 


-19.7 


1381.646 


0.646 


12.7 


-12.7 


1382.668 


0.738 


5.5 


-11.5 


1383.645 


0.826 


14.2 


-12.4 


1384.625 


0.914 


15.2 


-15.1 


1385.687 


0.010 


34.4 


-21.3 


1386.650 


0.096 


87.7 


-16.1 
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Table 5 — Continued 



Date 


Orbital 


V r 


W x 
( A ) 


( ±1 J L) — Z,4oU,UUU ) 


Phase 


(km s ) 


1388.583 


0.270 


54.8 


-11.8 


1390.664 


0.457 


40.5 


-19.0 


1391.602 


0.541 


11.3 


-18.7 


1392.619 


0.633 


10.8 


-13.2 


1393.638 


0.724 


11.7 


-13.1 


1394.615 


0.812 


17.2 


-13.0 


1421.765 


0.253 


71.7 


-12.1 


1425.784 


0.614 


19.5 


-14.3 


1426.745 


0.700 


22.9 


-12.4 


1427.748 


0.790 


41.3 


-12.4 


1428.718 


0.878 


36.7 


-13.1 


1429.713 


0.967 


37.7 


-20.9 



Table 6. Features Associated with the Massive Companion 



Date 
(HJD-2,450,000) 


Orbital 

Phase 


K(Si III) 
(km s- 1 ) 


y r (He II) 
(km s -1 ) 


iy A (He II) 

(A) 


1373.693 


0.931 




51.3 


-0.93 


1373.758 


0.937 




51.1 


-0.92 


1374.655 


0.018 




67.2 


-0.99 


1375.709 


0.113 




63.3 


-0.95 


1379.718 


0.473 




80.2 


-0.63 


1381.646 


0.646 




51.1 


-0.77 


1382.668 


0.738 


-97.8 


10.3 


-0.63 


1383.645 


0.826 


-59.7 


39.0 


-0.74 


1384.625 


0.914 




71.8 


-0.74 


1385.687 


0.010 




51.6 


-0.95 


1386.650 


0.096 




51.4 


-0.96 


1388.583 


0.270 


3.6 


122.0 


-0.66 


1390.664 


0.457 




80.8 


-0.60 


1391.602 


0.541 




7.8 


-0.74 


1392.619 


0.633 




24.7 


-0.81 


1393.638 


0.724 


-84.6 


15.6 


-0.78 


1394.615 


0.812 


-59.7 


1.1 


-0.58 


3152.930 


0.871 




86.8 


-0.59 


3153.825 


0.951 




59.2 


-0.70 


3154.928 


0.050 




72.9 


-0.71 


3156.922 


0.229 


33.0 


106.5 


-0.29 


3157.932 


0.320 


73.9 


123.5 


-0.34 


3158.928 


0.410 




72.1 


-0.58 


3161.915 


0.678 




60.4 


-0.34 


3162.914 


0.768 


-82.6 


23.2 


-0.51 


3163.892 


0.856 




35.3 


-0.52 


3164.903 


0.947 




54.0 


-0.66 
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Table 7. Summary of Radial Velocity Fits 



Element 


K(Si III) 


V r (Re II A4686) 


\4(Ho! wings) 


K (km s- 1 ) 


59 ±5 


36 ±6 


30 ±3 


7 (km s -1 ) 


-20 ±5 


64 ±4 


49 ±2 


r.m.s. (km s -1 ) 


22 


21 


13 
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■1000 1000 

VELOCITY (km s 1 ) 



Fig. 1. — The orbital phase variations in the Si IV A4088 absorption line in the FEROS 
spectra of RY Scuti are shown in linear plots (top panel) and as a gray-scale image (lower 
panel). The intensity in the gray-scale image is assigned one of 16 gray levels based on its 
value between the minimum (dark) and maximum (bright) observed values. The intensity 
between observed spectra is calculated by a linear interpolation between the closest observed 
phases (shown by arrows along the right axis). The feature centered at km s _1 is Si IV 
A4088, and to the right is N III A4097, H5 A4101 (both absorption and nebular emission), 
and Si IV A4116 (far right edge). 
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Fig. 2. — The orbital phase variations in the He I A4387 singlet shown in the same format 
as Fig. [TJ This feature strengthens at conjunctions and develops a blueshifted absorption 
feature after <fi = 0.5. Traces of the radial velocity curve of the fainter massive companion 
are seen especially after <p = 0.5, however, this feature is not used in the radial velocity 
analysis due to the aforementioned blueshifted absorption feature. The fainter absorption 
feature just to the left is N III A4379. 
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Fig. 3. — The orbital phase variations in the He I A4471 triplet shown in the same format as 
Fig. [TJ This line has a stronger blueshifted absorption feature after <fi = 0.5 that transforms 
into a narrow and long-lived absorption line. There is also some evidence of redward emission 
at the conjunctions that hints at a P Cygni profile shape. 



-34- 



m 

DC 

o 



0.0 



LU 

tt> 0.5 
< 

X 
Q_ 



1.0 



0.0 



0.5 



1.0 




Si III X4552 



20% 





Vv 



I 'I' 



ft 



2000 



■1000 o 1000 

VELOCITY (km s 1 ) 



■s 



2000 



Fig. 4. — The orbital phase variations in the Si III A4552 line in the same format as Fig.[U All 
of the triplet Si III A4552, 4567, 4574 (at 0, 1004, and 1459 km s _1 , respectively) components 
show dramatic strengthening at the conjunctions and display a broad, shallow feature moving 
in anti-phase as expected for the massive companion. The absorption feature at —726 km s _1 
is He II A4541. 
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Fig. 5. — The orbital variations in the Ha profiles after subtraction of the nebular component 
(based upon the KPNO coude feed set of spectra). The dashed line at <fi = 0.15 is a single 
HST spectrum of the central binary exclusive of the surrounding nebula, and the good match 
to our spectra at that phase indicates that the nebular subtraction method is reliable. The 
left hand panel shows the profiles normalized to the observed continuum (lower near = 0.0 
and = 0.5) while the right hand panel shows the emission rescaled to a constant reference 
continuum. The net Ha profile is approximately stationary and has high velocity wings that 
move in anti-phase to the supergiant's radial velocity curve. 
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Fig. 6. — The orbital phase variations of the broad and weak He II A4686 emission line. It 
appears brighter at <fi = 0.0 and shows a hint of the anti-phase motion associated with the 
massive companion. The nebular emission feature at 1000 km s -1 is [Fe III] A4702. 



0.0 



0.5 

ORBITAL PHASE 



1.0 



Fig. 7. — The radial velocity curve and measurements for several spectral components (where 
= 0.0 corresponds to supergiant superior conjunction). The filled circles represent radial 
velocities of the supergiant (with typical errors of 10 km s -1 ) while the large amplitude solid 
line shows the derived velocity curve for the circular solution and the elliptical solution is 
shown by the dotted curve (Table HJ). The open circles illustrate the radial velocities for 
Si III A4552 (with typical errors of 30 km s _1 ) and the smaller amplitude solid line shows 
the constrained fit (Table [7]). The diamonds represent the radial velocities for the broad 
He II A4686 emission (with typical errors of 50 km s" 1 ) while the plus signs indicate radial 
velocities for the Ha wings (with typical errors of 2 km s _1 ). 
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Fig. 8. — Tomographically reconstructed spectra for the supergiant and the environment of 
the massive companion in the spectral range 3950 — 4250 A plott ed using normalized flux. 
These are compared with two single star spectra from the atlas of IValdes et al.l (120041 ) . The 
spectra are offset in flux for clarity. The reconstructed spectra display some "rippling" near 
some strong features due to residual emission with non-orbital motions. 
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Fig. 9. — Tomographically reconstructed spectra in the range of 4450 — 4750 A (in the same 
format as Fig. [S]) . 
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Fig. 10. — Cartoon model of the system as seen from above the orbital plane with the ob- 
server's orbital phases noted on the periphery and a scale bar denoting lOR® in the lower 
right of the diagram. The Roche-filling supergiant appears on the left while the massive com- 
panion is shown on the right surrounded by a dense accretion disk with the outer boundary 
shown as the dotted line. The mass loss from the supergiant primarily occurs in two regions: 
the LI region between the stars and the L2 region on the left side of the supergiant. These 
are delineated by the dashed lines. The thick solid line represents the classical Roche lobe 
overflow stream from LI to the disk. The tick mark on the axis joining the centers of the 
two stars marks the center of mass, and the arrows at the centers of both stars indicate their 
orbital velocities. The LI and L2 points lie along the axis joining the two stars at the points 
indicated. The L3 point lies off the plot to the right (62i? from the center of the massive 
companion). 
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Fig. 11. — A representation of the apparent increases in equivalent width due to the changing 
continuum lev els in the eclipsin g binar y. The small plus signs mark the inverse fluxes in the 
V^-band from iDjurasevic et all (120011 ). and these show the predicted trend in equivalent 
width for any uneclipsed emission source. The filled circles show the normalized variation 
in Ha equivalent width that appears to follow the predicted curve (indicating that the Ha 
source is not significantly eclipsed). The open circles show the same for He II A4686. In 
this case, the absence of brightening at <p = 0.5 suggests that the He II emission is partially 
occulted/eclipsed then. Conservative estimates of the relative equivalent width errors for 
Ha are 4% and those for He II A4686 are between 4% and 10%. 



